Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 20 October 2009 (MN WI^ style file v2.2) 



Two populations of metal-free stars in the early Universe 



Thomas H. Greif" and Volker Bromm 

Department of Astronomy, University of Texas, Austin, TX 78712, USA 



20 October 2009 



o 
o 

(N 

-(— > 
o 

O 

o 

> 

o 
o 

6 



% 



ABSTRACT 

We construct star formation histories at redshifts 2^5 for two physically distinct 
populations of primordial, metal-free stars, motivated by theoretical and observational 
arguments that have hinted towards the existence of an intermediate stellar generation 
between Population III and Population I/II. Taking into account the cosmological 
parameters as recently revised by the Wilkinson Microwave Anisotropy Probe after 
three years of operation, we determine self-consistent reionization histories and discuss 
the resulting chemical enrichment from these early stellar generations. We find that 
the bulk of ionizing photons and heavy elements produced at high redshifts must have 
originated in Population II. 5 stars, which formed out of primordial gas in haloes with 
virial temperatures > 10* K, and had typical masses > 10 M0. Classical Population III 
stars, formed in minihaloes and having masses > 100 Mq, on the other hand, had only 
a minor impact on reionization and early metal enrichment. Specifically, we conclude 
that only ~ 10 per cent by mass of metal-free star formation went into Population III. 
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1 INTRODUCTION 

Understanding the formation and properties of the first 
stars, formed out of pure H/He gas at the end of the cosmic 
dark ages, is one of the key problems in modern cosmol- 

, ogy (e.g. Barkana & Loeb 2001; Bromm & Larson 2004). 

' Within the standard A cold dark matter (ACDM) model of 
structure formation (e.g. Spergel et al. 2006), the first stars 

, are predicted to have formed in 'minihaloes' of total mass 
> 10*^ Mq at redshifts ^ > 20 (Haiman, Thoul & Loeb 1996; 
Tegmark et al. 1997; Yoshida et al. 2003). Their emergence 
signalled a fundamental transition in cosmic history, from 
an initial state of simplicity to one of increasing complexity 
(e.g. Ciardi & Ferrara 2005). Specifically, the first stars en- 
riched the intergalactic medium (IGM) with the first heavy 
elements (Gnedin & Ostriker 1997; Madau, Ferrara & Rees 
2001; Mackey, Bromm & Hernquist 2003; Furlanetto & Loeb 
2003, 2005; Scannapieco, Schneider & Ferrara 2003; Bromm, 
Yoshida & Hernquist 2003; Yoshida, Bromm & Hernquist 
2004; Schneider et al. 2006). Secondly, they contributed to 
the reionization of the Universe, as inferred from the optical 
depth to Thomson scattering measured by the Wilkinson 
Microwave Amsotropy Probe ( WMAP; Spergel et al. 2006) , 
due to their efficient production of ionizing photons. To elu- 
cidate these feedback effects, we have to ask: What were the 
properties of the first generations of metal-free stars? 
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Numerical simulations imply that star formation in- 
side 'minihaloes' relied on molecular hydrogen cooling, 
and resulted in typically very massive so-called Popula- 
tion III (Pop HI) stars with characteristic masses > 100 M© 
(Bromm, Coppi & Larson 1999, 2002; Abel, Bryan & Nor- 
man 2000, 2002; Nakamura & Umemura 2001). Minihaloes 
are characterized by virial temperatures below 10* K, and 
the star forming gas did not go through an ionized phase, 
prior to the late stages of the protostellar accretion process. 
In environments, however, where such ionization did occur, 
the thermal evolution of the primordial gas is predicted 
to have been markedly difi'erent (see Johnson & Bromm 
2006a,b for details) . This difference arises because the hydro- 
gen deuteride (HD) molecule provides an additional cooling 
channel, allowing the metal-free gas to cool to temperatures 
below those possible solely via H2 cooling (e.g. Flower et al. 
2000; GaUi & Palla 2002; Nakamura & Umemura 2002). Gas 
fragmentation in cases where HD cooling is important could 
therefore have resulted in metal-free stars that were some- 
what less massive than 'classical' Pop HI stars. These stars, 
having typical masses > 10 M0 , have been termed Popu- 
lation II. 5 (Pop II. 5), indicating their transitional status in 
terms of typical mass and formation epoch between Pop HI 
and Pop I/II (Mackey et al. 2003; Johnson & Bromm 2006a). 
Connected to this, Tumlinson, Venkatesan & ShuU (2004) 
and Tumlinson (2006a, b), together with other groups, have 
argued that the observed chemical abundances in metal-poor 
Galactic halo stars are best explained by assuming that the 



© 0000 RAS 



2 T. H. Greif and V. Bromm 



first stars were somewhat less massive, different from tlie 
theoretical expectation (see Section 4). Our model could 
provide a natural explanation for such an intermediate mass 
scale, both in terms of the physics of primordial gas cooling 
and in terms of the cosmological environment where such 
stars could have formed. 

The conditions necessary for HD cooling to become im- 
portant were possibly met in a number of high- 2; environ- 
ments (Johnson & Bromm 2006a). Among them arc fossil 
H II regions, left behind after the death of a Pop III star (e.g. 
Nagakura & Omukai 2005; Susa & Umemura 2006; Johnson 
& Bromm 2006a,b; Abel, Wise & Bryan 2006), and shocks 
induced either by the first supernova (SN) explosions (e.g. 
Machida et al. 2005), or as a result of hierarchical struc- 
ture formation during the assembly of the first dwarf galax- 
ies (e.g. Uehara & Inutsuka 2000; Nakamura & Umemura 
2002). In this paper, we consider in particular the latter 
possibility, assuming that Pop II. 5 stars formed during the 
virialization of primordial haloes with virial temperatures 
> 10* K, which is sufficient to ionize the collapsing gas. 

In the following, we construct separate cosmic star for- 
mation histories for the two distinct populations of metal- 
free stars (Section 2). We then compare the ionizing photon 
and heavy element production of our model to observations 
(Section 3), and conclude by summarizing our results and 
discussing the relevance of our work in light of recent devel- 
opments (Section 4). 

Throughout this paper, we adopt a ACDM model of 
hierarchical structure formation, with cosmological param- 
eters consistent with the WMAP 3-yr results (Spergel et 
al. 2006), i.e. density in matter fim = 1 — ^a = 0.24, 
density in baryons Sib = 0.042, Hubble parameter h — 
Ho/ (100 km s"^ Mpc"^) = 0.73, tilt n^ = 0.95 and a nor- 
malization of the power spectrum to ag — 0.74. For calcu- 
lations involving the latter, we use the transfer function of 
Eisenstein & Hu (1999). 



2 COSMIC STAR FORMATION HISTORY 

2.1 Population III 

As the very first generation of stars to form out of primor- 
dial gas. Pop III has been investigated extensively over the 
last decade (e.g. Barkana & Loeb 2001; Bromm & Larson 
2004; Ciardi & Ferrara 2005; Glover 2005). Already prior to 
the first numerical simulations of early structure formation, 
molecular hydrogen was identified as the primary coolant 
for pure H/He gas (e.g. Peebles & Dicke 1968; Couchman & 
Rees 1986; Tegmark et al. 1997). This view has been con- 
firmed by detailed simulations since 1996 (e.g. Haiman et al. 
1996; Abel et al. 1998, 2000; Bromm et al. 1999, 2002; Fuller 
& Couchman 2000; Nakamura & Umemura 2001; Yoshida 
et al. 2003, 2006). Within variants of the CDM model, the 
first stars are predicted to have formed in minihaloes of to- 
tal mass ~ 10^ M0 at redshifts z > 20. Such minihaloes 
have virial temperatures in excess of a minimum value, 
Tmin — 10"^ K, required for collapsing gas to cool in less than 
a Hubble time (e.g. Haiman et al. 1996; Tegmark et al. 1997, 
Yoshida et al. 2003) . The evolution of gas inside minihaloes, 
driven by H2 cooling, leads to the formation of protostel- 
lar fragments with typical masses > 100 M© (Bromm et al. 



1999, 2002; Abel et al. 2000, 2002; Yoshida et al. 2006). We 
adopt this characteristic mass scale, and assume throughout 
this paper that Pop HI has an extremely top-heavy initial 
mass function (IMF) with cutoffs at 100 and 500 M© , and a 
standard Salpeter (1955) slope. Consistent with simulation 
results (Abel et al. 2000, 2002; Bromm et al. 2002) and re- 
cent semianalytic treatments predicting a low Pop HI star 
formation efficiency (e.g. Haiman & Bryan 2006; Choudhury 
& Ferrara 2006), we assume /, ~ 0.001, which corresponds 
to > 1 stars per halo. However, since uncertainties remain, 
we will reconsider our choice of this crucial quantity in later 
sections. 



2.1.1 Population III Star Formation 

A straigthforward semianalytic approach to determine the 
star formation rate (SFR) for Pop HI is given by the Press- 
Schechter (PS) formalism (Press & Schechter 1974), provid- 
ing an expression for the abundance of dark matter haloes. 
An improved fit to simulation results is obtained by assum- 
ing ellipsoidal instead of spherical collapse, introduced by 
Sheth, Mo & Tormen (2001). In quantitative studies (e.g. 
Heitmann et al. 2006; Reed et al. 2006), the Sheth- Tormen 
(ST) approach has proven to be significantly more accurate, 
with maximum deviations of ~ 50 per cent, leading us to 
adopt this formalism for all further calculations. By com- 
parison, the PS mass function underpredicts the abundance 
of haloes at very high redshifts by almost an order of mag- 
nitude, although our qualitative results remain unaffected. 
Using the ST approximation, we can estimate the num- 
ber of dark matter haloes, N{M,z), per unit mass and co- 
moving volume at a given redshift: 



= nsT{M,z) 



(1) 



dN{M,z) 

dMdV 

Since only a small fraction of haloes is massive enough to 
fulfil the cooling condition imposed by Tvir > Tmin = 10^ K, 
we must adjust the limits of integration appropriately when 
calculating the collapsed fraction of mass, Fcoi{z), available 
for Pop HI star formation: 



fc 



Pm 



dMMnsT{M,z) 



(2) 



where pm is the total mass density of the background uni- 
verse. The limits in the integral, Mmin ~ Afcrit(T'min) and 
Afmax = Afcrit(rmax), are related to the virial temperature 
by: 



.,_,„. M.(iL)-""(Jj^)-(i±i 



-3/2 



(3) 



where p is the mean molecular weight, i.e. 1.2 for neutral 
and 0.6 for ionized, primordial gas. Since the Universe is 
mostly neutral at high redshifts, we choose p — 1.2 for all 
further calculations. The lower mass limit in the integral cor- 
responds to the minimum virial temperature below which 
haloes cannot cool efficiently. We set the upper limit to 
Tmax = 10'' K, as haloes with masses > 10* Mq go through 
a more complicated cooling process, leading to the forma- 
tion of stars with characteristic masses j^ 10 Mq, termed 
Pop II. 5. This novel generation of stars will be discussed in 
Section 2.2. 
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Applying the above criteria, we obtain a preliminary 
form for the Pop III SFR: 






dFcoi(z) 



dz 



(4) 



The solid line in Panel (b) of Fig. 1 shows the result, imply- 
ing that Pop III star formation peaks around z ~ 10, and 
declines sharply towards z ~ 5, as an increasing number 
of minihaloes merge into larger haloes with masses above 
Afmax- We emphasize that this is not the final SFR, as we 
have neglected all forms of feedback. In the next few sec- 
tions, we will therefore investigate the impact of relevant 
feedback mechanisms on Pop III. 



2.1.2 Feedback Mechanisms 

At some point during primordial star formation, negative 
feedback effects must have become important, eventually 
terminating the Pop III mode, as stars in the present-day 
Universe are characterized by much higher metallicities and 
typically have low masses. Perhaps the most relevant neg- 
ative feedback mechanism is heating due to photoioniza- 
tion (photoheating), since massive Pop III stars produce 
vast amounts of ionizing photons (e.g. Tumlinson & ShuU 
2000; Bromm, Kudritzki & Loeb 2001b; Schaerer 2002), 
with escape fractions around unity (Whalen, Abel & Nor- 
man 2004; Kitayama & Yoshida 2005; Alvarez, Bromm & 
Shapiro 2006a). Related to this, other radiative feedback 
mechanisms have also been considered (for a recent review, 
see Ciardi & Ferrara 2005). 

Of particular importance is the impact of radiation 
in the Lyman-Werner (LW) bands on neighbouring mini- 
haloes, which might readily dissociate Kg and suppress ef- 
ficient cooling (e.g. Haiman, Rees & Loeb 1997; Omukai & 
Nishi 1999; Ciardi, Ferrara & Abel 2000; Haiman, Abel & 
Rees 2000; Machacek, Bryan & Abel 2001; Ricotti, Gnedin 
& ShuU 2001, 2002a,b; Kitayama et al. 2001; Glover & 
Brand 2001; Yoshida et al. 2003; Omukai & Yoshii 2003; 
Mesinger, Bryan & Haiman 2006; Maclntyre, Santoro & 
Thomas 2006). Despite the extensive effort, no consensus 
has yet been reached on the overall significance of this ef- 
fect. Alternatively, positive feedback due to an early X-ray 
background produced by relic black holes has been suggested 
(e.g. Haiman et al. 2000; Machacek, Bryan & Abel 2003; 
Glover & Brand 2003). Pop HI stars in the range 100 - 140 
and 260 - 500 Mq are expected to collapse directly to black 
holes after their main sequence lifetimes of < 3 Myr (Heger 
et al. 2003), emitting X-rays during the subsequent accre- 
tion of gas. Other possible sources include X-ray binaries 
and SN remnants (e.g. Natarajan & Almaini 2000; Helfand 
& Moran 2001). The buildup of an early X-ray background 
could enhance the production of free electrons in neighbour- 
ing minihaloes, thereby increasing the H2 fraction and con- 
sequently the cooling rate. This would result in an increase 
of Pop HI star formation and could possibly compensate for 
negative feedback caused by photoheating or LW radiation. 

Based on the revised value of the optical depth to 
Thomson scattering (e.g. Spergel et al. 2006), we find that 
Pop HI must have been terminated fairly rapidly (see Sec- 
tion 3). This leads us to concentrate on negative feedback, 
and in the following we will discuss the impact of the LW 



background, and subsequently of photoheating, on primor- 
dial star formation. 



2.1.3 Lyman-Werner Feedback 

The high masses and surface temperatures of Pop HI stars 
give rise to yields of the order of 10^ ionizing photons per 
stellar baryon (e.g Bromm et al. 2001b; Schaerer 2002). Nat- 
urally, the flux in the LW bands is also high, dissociating H2 
in neighbouring haloes and reducing the fraction of gas that 
can cool efficiently. To describe this feedback mechanism, 
we must find a connection between the background LW flux 
and the cooling rate of a halo, and also between the mass in 
stars and the associated photon yields. 

For simplicity, we assume that all emitted LW pho- 
tons escape into the IGM, and are thus available for H2- 
dissociation in a nearby minihalo, but account for the red- 
shifting of photons out of the LW bands. This approach 
clearly overestimates the LW background, yet will turn out 
to be instructive when discussing the results. The buildup of 
a homogeneous LW background proceeds over a very short 
period of time, since the comoving mean free path of a LW 
photon at z = 30 is ~ 10 Mpc (Mackey et al. 2003). Of- 
ten, the presence of such a pervading background has been 
assumed in investigating the cooling process of a minihalo 
under LW irradiation. Specifically, Machacek et al. (2001) 
have derived an analytic fit to their simulation results, giv- 
ing an estimate of the gas fraction that cools and reaches 
high densities as a function of halo mass and LW flux: 



/Lw(Af,J2i)=0.061n 



M/Mg 



1.25 X 105 + 8.7 X 105R°,4'' 



, (5) 



where F21 = 47rJ2i is the flux in the LW bands, such that: 
Jlw ~ J21 X 10^^^ erg s^^ cm^^ Hz^^ sr^^. This formula 
explicitly shows that H2 cooling rates rise towards higher 
halo masses, enabling them to withstand larger amounts of 
LW flux. Averaging this fraction over the halo mass range 
hosting Pop HI, and 'normalizing' it to the value obtained 
for an absent LW background, yields the mass fraction of 
gas available for star formation as a function of z and J21 : 



jM^x dMMnsT(M, 2)/lw(M, J21) 

LW(J21,Z) = — p ; 

jM„.x dMMnsT(M, 2)/lw(M, J21 = 0) 



,(6) 



where the normalisation accounts for the fact that Machacek 
et al. (2001) imply a definition for the star formation effi- 
ciency different from our choice of a constant /, , indepen- 
dent of halo mass. Thus, the SFR at any given redshift is 
reduced relative to the no-feedback SFR of equation (4) by 
precisely this factor: 



*UJ21,z)=/lw(J21,z)*.(z) 



(7) 



where J21 is implicitly a function of redshift as it is linked 
to the mass already incorporated into stars. To complete 
the feedback cycle, we must therefore connect the LW flux 
with the comoving density in stars, p*{z). For simplicity, we 
assume that the ratio of mean LW frequency to the relevant 
bandwidth is of the order of unity: 



Jhwiz) 



he / \ /, n3 

nLWP*{z) (1 + Z) 



47rmH 



(8) 
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where r/LW is the number of photons emitted in the LW 
bands per steUar baryon, and: 

di 



*,(2 



dz' 



dz' 



(9) 



such that Az = Zr — z corresponds to the maximum dis- 
tance a LW photon can travel before it is redshifted out 
of the LW bands. Since all stellar populations contribute 
to the LW background, possibly with distinct values for 
J7LW, the star formation histories of Pop IL5 and Pop I/II 
are directly linked to Pop III. Schaerer (2002) has pro- 
vided time-integrated spectral properties of metal-free stars 
over a wide range of masses, and our calculations of t^lw 
are based on this work. For Pop III, we integrate over the 
Salpeter- weighted mass range 100 - 500 M©, which results 
in 77LW — 10'*. As will be discussed in more detail later in this 
paper. Pop II. 5 is assumed to have masses > 10 M0, lead- 
ing us to adopt mass cut-offs at 10 and 100 M©, and again 
a standard Salpeter slope. Using the spectra of Schaerer 
(2002), we find r/LW ^ 2 x lO''. For Pop I/II, we adopt 
?7lw = 4 X lO"', which is comparable to the ionizing flux. The 
explicit SFRs used for Pop II. 5 and Pop I/II correspond to 
the feedback parameter KJ = 2/3, which will be discussed 
later. We note, however, that varying this parameter hardly 
affects our results. 

Applying the above treatment allows us to describe 
Pop III star formation in an evolving LW background. Our 
calculation is shown in Fig. 1 (dotted line), compared to 
the no- feedback case (solid line) . Somewhat surprisingly, the 
overall amplitude of the effect is not too high, although we 
have conservatively assumed that no LW photons are lost. 
On the other hand, we have not included the effect of bi- 
asing since the build-up of a homogenuous LW background 
is very rapid. The shape of the modified SFR is also quite 
interesting, and can be explained by a decomposition into 
two competing processes: As the LW background increases, 
the Pop III SFR tends to decline, yet as haloes become more 
massive, they can sustain higher LW fluxes, enabling ongo- 
ing Pop III star formation towards lower redshifts. To sum- 
marize, LW feedback might not be as efficient as previously 
assumed (e.g. Haiman et al. 1997; Mackey et al. 2003), yet 
for definitive conclusions one must apply detailed numerical 
simulations that self-consistently solve the hydrodynamics 
and radiative transfer. 



2.1.4 Photoheating 

The emergence of the first H 11 regions around star-forming 
haloes has an important effect on high-redshift star for- 
mation. The temperature of ionized bubbles is raised to 
~ 10* K, quenching star formation in affected regions, as 
primordial gas can no longer condense into the shallow grav- 
itational wells of minihaloes. Since Pop III stars produce co- 
pious amounts of ionizing photons, with escape fractions of 
the order of unity (e.g. Alvarez et al. 2006a), and haloes are 
strongly biased at these high redshifts, photoheating could 
in fact be a dominant feedback mechanism for Pop III (e.g. 
Kramer & Haiman 2006). To quantify this feedback, we be- 
gin with given ionization histories, and assume that Pop III 
star formation is suppressed by a factor equal to the volume 
filling fraction of ionized regions, Qion(z), weighted by the 
mean correlation between haloes. With this treatment, we 



1.000 r 




Figure 1. (a) Buildup of the LW background in units of 
10^^^ erg s^^ cm~^ Hz~^ sr"'^. We here consider the LW radia- 
tion produced by all stellar populations, (b) The SFR for Pop III, 
neglecting all feedback effects (solid line), and after inclusion of 
LW feedback (dotted line). 



might overestimate the strength of the feedback, as already- 
collapsed haloes are subject to photoevaporation rather than 
infall suppression due to photoheating, the former being a 
much slower process (e.g. Shapiro, Iliev & Raga 2004; Iliev, 
Shapiro & Raga 2005). However, since the collapse of neigh- 
boring minihaloes is synchronized only to within ~ lO'^ yr, 
compared to propagation times of the first Hll regions of 
~ 10^ yr (e.g. Alvarez et al. 2006a), we expect that most 
haloes have not yet reached high enough densities to be sub- 
ject to photoevaporation. 

To explicitly determine the effect of biasing, we use the 
well-known linear correlation function 5hh(Af , R, z) between 
two haloes of mass M at a comoving distance R: 



ii,i,{M,R,z) = h'{M,z)£,^m{R) , 



(10) 



where 6(M, z) is the linear bias of a halo of mass M and 
Cmm(i?) is the mass correlation function (e.g. Mo & White 
2002). Since haloes are distributed according to the ST func- 
tion, we can determine the average bias at any given redshift: 

jM„„ ^MMh{M, z)nsT{M, z) 



h{z) = 



J^'"-^'^dMMnsT{M,z) 



(11) 



which replaces b{M,z) in equation (10). The second step 
consists of linking the volume filling fraction of ionized re- 
gions to the typical size of an ionized bubble, which we do 
by assuming that the volume of a typical H 11 region is equal 
to Qion(z) divided by the current number density of haloes, 
n{z): 



^nR\z) 



n{z) 



(12) 



With these preparations, we can ascertain the probabil- 
ity pion(z) that a freshly collapsed halo lies within an already 
existing Hll region: 
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Pion{z) 



,{Z)[1 + ^M.{Z)] 



(13) 



where the halo correlation function (,hh{z) adopts the inter- 
pretation of an excess probability. This enables us to write 
the modified SFR with respect to the no-feedback SFR as: 



*:(z) = [l-pion(2)]*,(z) 



(14) 



To calculate the SFRs we must supply reionization histo- 
ries, and choose Gaussian functions of the form Qion{z) = 
g-(2-6) /w ^ where w = 57, 25 and 5 corresponds to early, in- 
termediate and late reionization, respectively. For all three 
cases, reionization ends at z = 6, thus being in agree- 
ment with observations suggesting that the Universe was 
completely ionized by z>6, e.g. the absence of the Gunn 
& Peterson (1965) trough in the spectra of high-redshift 
quasars, and considerations of the thermal history of the 
IGM (e.g. Theuns et al. 2002; Hui & Haiman 2003). Ad- 
ditionally, the three reionization histories reproduce optical 
depths of r = 0.12, 0.09 and 0.06, thus spanning the cur- 
rent error bars of the WMAP measurement. For the three 
different values of w as mentioned above. Panel (a) in Fig. 2 
shows early (dotted line), intermediate (dashed line) and 
late (dot-dashed line) reionization, whereas Panel (b) traces 
the modified SFRs derived from these reionization histories. 
As is evident from the plots, the difi'erent epochs of reioniza- 
tion lead to vastly different SFRs, with efficiencies varying 
by more than one order of magnitude. This is primarily an 
effect of halo biasing, amplifying the impact of small vol- 
ume filling fractions Q\on[z). We tentatively conclude that 
Pop III ends very rapidly once a small number of minihaloes 
have formed stars. 

At this point, it is important to note that we have not 
yet derived self-consistent reionization and star formation 
histories. However, it is instructive to outline the range of 
possibilities and defer conclusions about the most likely and 
self-consistent scenario to later sections. 



2.2 Population II. 5 

The term Pop III has traditionally been used to denote stars 
forming in extremely low metallicity environments, indepen- 
dent of the associated halo mass. In the following, we will 
deviate from this terminology, and instead distinguish be- 
tween two populations of metal-free stars according to the 
environments in which they form. We thereby propose an in- 
termediate stellar generation between Pop III and Pop 1/11 
(see Mackey et al. 2003; Johnson & Bromm 2006a). 

As discussed in Section 2.1, haloes in the range Tvir — 
10"^ - 10* K cool almost exclusively by molecular hydrogen 
transitions, and produce stars with an extremely top-heavy 
IMF. Yet gas assembled in more massive haloes goes through 
an ionized phase and experiences a different cooling process 
(Johnson & Bromm 2006a, b). In this process, HD cooling 
allows the gas temperature to reach that of the cosmic mi- 
crowave background (CMB) within a Hubble time. By reach- 
ing the CMB limit, these haloes produce stars with typical 
masses > 10 M0 , about an order of magnitude less than is 
the case for Pop 111 (Johnson & Bromm 2006a). Associating 
this intermediate generation of stars with metal-free haloes 
above Tvir = 10* K allows us to derive a SFR for Pop 11.5. 



1.0 

0.8 

c 0.6 

0.4 

0.2 

0.0 

10"^ 

IQ-* 



" IQ- 

Li_ 
(/) 

10" 



: 


\ 
\ \ 




(0): 


7 


\ \ 
\ \ 

\ \ 
\ 




- 




\ 




- 










/ 


/ 


\ 


ib): 


; / 




t \ 





10 



15 20 25 30 

z 



Figure 2. (a) Early (dotted line), intermediate (dashed line) and 
late (dot-dashed line) reionization, here applied to calculate the 
suppression of Pop III star formation by photoheating. (b) Com- 
parison of the modified SFRs to the no- feedback case (solid line). 



2.2.1 Population II. 5 Star Formation 

To obtain a preliminary SFR for Pop II. 5, we once again 
apply the ST formalism, yet adjust the limits of integration 
when determining the collapsed mass fraction by setting the 
lower limit to Mmin = Merit (Tmin = 10* Kj and the upper 
limit to infinity. This satisfies the conditions for their forma- 
tion sites, as these haloes go through an ionized phase in the 
course of their virialization (e.g. Johnson & Bromm 2006a). 
The SFR is computed according to equation (4), except that 
we use a different star formation efficiency, /« ~ 0.01. This 
choice refiects the intermediate nature of Pop 11.5 between 
Pop 111 and Pop 1/11, and will be further justified in later 
sections. The solid line in Panel (b) of Fig. 3 shows the re- 
sulting SFR, implying that Pop 11.5 produces stars about an 
order of magnitude more efficiently than Pop HI (see Fig. 2), 
but also peaks at later times. 



2.2.2 Feedback Mechanisms 

Due to the deep potential wells of typical Pop 11.5 haloes 
(>10* Mq), and the fact that cooling by atomic hydro- 
gen sets in, photoheating is no longer a dominant feedback 
mechanism. Specifically, Dijkstra et al. (2004) have shown 
that suppression of gas infall is not as efficient in the high- 
redshift Universe as it is for z<6. We performed calcula- 
tions including infall suppression in ionized regions of the 
IGM and verified this conclusion. Efficient self-shielding of 
ionized gas assembled in more massive haloes also renders 
LW feedback unimportant (see Johnson & Bromm 2006b). 
Consequently, Pop 11.5 can only be terminated by chemi- 
cal feedback, as multiple SNe in a host galaxy combine to 
form so-called 'superwinds', expanding into the IGM and 
enriching neighbouring haloes with metals (e.g. Madau et 
al. 2001; Mori, Ferrara & Madau 2002; Furlanetto & Loeb 
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2003, 2005; Wada & Venkatesan 2003). In the next section, 
we will therefore introduce a semianalytic model describing 
chemical enrichment, which subsequently leads to the ter- 
mination of Pop II. 5. 

2.2.3 Chemical Feedback 

Once the metallicity of the IGM exceeds a critical level, the 
formation of primordial, metal-free stars is suppressed and 
star formation switches to a more conventional mode (e.g. 
Omukai 2000; Bromm et al. 2001a; Schneider et al. 2002). 
The precise value of this 'critical metallicity' is still heav- 
ily debated, as some authors argue that dust production in 
the early Universe might be responsible for lowering Zcrit 
to almost 10~^ Zq (e.g. Schneider et al. 2006, and refer- 
ences therein). However, due to remaining uncertainties we 
choose the more robust value Zcth = 10~^'^ Zq, which is 
based on the physics of fine structure cooling, and does not 
depend on whether dust is present or not, and if so, what 
its composition might be (e.g. Bromm et al. 2001a; Bromm 
& Loeb 2003; Santoro & ShuU 2006). Pop III stars in the 
range 140 - 260 Mq explode as Pair-Instability Supernovae 
(PISNe), with substantial metal yields y ~ 0.5 (Heger et 
al. 2003), yet due to their rare nature they are likely ineffi- 
cient polluters (also see Section 3.3). We therefore suggest 
that Pop II. 5 is largely responsible for enriching the IGM 
above the critical threshold and quenching metal-free star 
formation. 

In modelling the propagation and distribution of metals, 
we closely follow the semianalytic treatment of SNe winds 
performed by Furlanetto & Loeb (2005). Utilizing the well- 
known Sedov (1959) solution of an expanding bubble in a 
uniform background, they find an expression for the ratio of 
enriched mass to (total) halo mass: 

-3/5 

ri(M) = 27 K^ ft' "^Eli"^ 



M 



1010 Mq 



100 M© 

2/5 



(^) 



-3/5 



(15) 



where K^ is a free parameter describing radiative losses, £51 
is the typical energy emitted by a SN in units of 10'^ erg and 
iOs-n is the mass in stars per SN (~ 20 Mq for our assumed 
Pop II. 5 IMF). By integrating r]{M) over the ST function 
to determine the fraction of space with metals, and weight- 
ing it with the mean bias at a given reshift, Furlanetto & 
Loeb derive a probability function Ppris(z) for freshly col- 
lapsed gas to be pristine. We identify this quantity with the 
fraction of gas that has just collapsed and is available for 
Pop II. 5 star formation. In doing so, we implicitly assume 
that the mean metallicity of the fraction of space reached by 
an expanding superwind lies well above Zc-At. We will justify 
this assumption below. 

The properties of these SNe follow directly from our 
assumed IMF. Stars in the range ~ 10-25 Mq explode 
as conventional core-collapse SNe with energies i?5i ~ 1, 
and metal yields y ~ 0.05, while stars with masses above 
~ 40 Mq collapse directly to black holes (Heger et al. 2003) . 
The fate of stars with masses 25 - 40 Mq seems to be not 
so unambiguous, as those with low core angular momenta 
end their lives as faint SNe, while those with higher core 
angular momenta may explode as hypernovae (HNe), with 



energy outputs E^i > 10 (Heger et al. 2003). As a reasonable 
compromise, we therefore assume that all stars in the range 
25 - 40 Mq explode as HNe with typical energies Er,i = 
10, although variations herein do not significantly affect our 
results. Averaging the energy output and metal yield over 
the IMF results in S51 ~ 2 and y ~ 0.04. 

With this information, we can determine the metal con- 
tent of a region affected by a typical superwind. Explicitly, 
the metallicity of an enriched region is given by the ratio of 
mass in metals M^ to the mass of the expanding bubble: 

The mass in metals can be calculated from the total mass 
of the star-forming region by means of the star formation 
efficiency and the metal yield, resulting in: 



Z = 



y 



2/5 p-3/5 



27K, 



M 
10i« Mq 



100 Mq 



3/5 



2/5 



10 



3/5 



(17) 



For a typical Pop II. 5 halo of 10 Mq, forming at z — 15, 
and with K^' set to the fiducial value 2/3, we find Z ~ 
3 X 10~^ Zq. This lies well above the critical metallicity and 
proves that our initial assumption was justified. 

To obtain a SFR that accounts for enrichment by su- 
perwinds, we simply multiply the original SFR with ppris(z), 
thereby selecting the pristine fraction of freshly collapsed 



gas: 



*1(2) =Ppris(2)**(z) 



(18) 



Using this semianalytic model, we compute star formation 
histories for three different radiation loss parameters, i.e. 
KvJ — 1, 2/3 and 1/2, corresponding to strong, intermedi- 
ate and weak chemical feedback, respectively. Panel (a) in 
Fig. 3 shows Ppris(z) for K„ = 1 (dotted line), K^J — 2/3 
(dashed line) and K^J — 1/2 (dot-dashed line), whereas 
Panel (b) shows the implied SFRs. As expected, higher ra- 
diative losses result in less chemical enrichment and more ef- 
ficient star formation. Also, the termination process is much 
slower than is the case for Pop III, due to the nature of 
the feedback: Wind hosts are heavily biased at these high 
redshifts, implying that enriching superwinds overlap and 
do not enrich large volume filling fractions, whereas the ef- 
fect of biasing only marginally inhibits the propagation of 
ionization fronts. 

2.3 Population I/II 

In order to compare the SFRs we have derived for Pop HI 
and Pop II. 5 to recent stellar generations, and also deter- 
mine the reionization history of the Universe, we will briefiy 
discuss our model for Pop I/II. Closely following the treat- 
ment of Pop II. 5, and maintaining consistency in terms of 
our ST approach, we assume that haloes above Tvir = 10* K 
form Pop I/II stars once they have been polluted by enrich- 
ing superwinds: 

*:(Z) = [1-Ppris(2)]*,(2), (19) 

where '^^{z) is the no- feedback SFR introduced in Sec- 
tion 2.2.1. Fig. 4 shows the resulting SFR for Pop I/II using 
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Figure 3. (a) Evolution oippy-iai^) for strong (dotted line), inter- 
mediate (dashed line) and weak (dot-dashed line) chemical feed- 
back, (b) Comparison of the no- feedback SFR for Pop II. 5 (solid 
line) to the SFRs modulated by Pp,-is{z). 



Li_ 




10 15 20 25 30 



Figure 4. Pop I/II star formation for intermediate chemical feed- 
back (dashed line), compared to an analytic fit by Springcl & 
Hernquist (solid line). 



K^x ~ 2/3, compared to an accurate fit of recent star forma- 
tion iiistory by Springel & Hernquist (2003) of tiie form: 



*.(«) ^K 



a exp [b {z - 



a — b + b exp [a (z 



0] 



(20) 



wliere K = 0.15 Mq yr"^ Mpc"^ a = 14/15, b = 3/5 
and Zm ~ 5.4. Note that this expression was derived for 
the WMAP first year parameters, and tiie true SFR would 
be somewhat lower. In the redshift regime of interest, i.e. 
z < 15, both plots are in good agreement, since at higher 
redshifts Pop II. 5 becomes the dominant mode. 

In the next section, we will evaluate the SFRs we have 
derived for Pop III, Pop II. 5 and Pop I/II by comparing 
the ionizing photon and metal production to observational 
constraints. 



3 OBSERVATIONAL SIGNATURES 

To test the validity of our models, we must compare our 
SFRs to observational signatures. This turns out to be a 
formidable task, as the sites of early star formation are the 
most distant sources of stellar radiation in the Universe, 
and all empirical constraints to date have been indirect. Yet 
the situation is rapidly improving by utilizing existing tele- 
scopes such as Spitzer, and pushing them to their limits 
(e.g. Barkana & Loeb 2001; Bromm & Larson 2004; Ciardi 
& Ferrara 2005; Beers & Christlieb 2005). 

One of the most important clues on early star formation 
is the WMAP measurement of the optical depth to Thomson 
scattering, recently revised to be r ~ 0.09 ±0.03 (Spergel et 
al. 2006). By providing an integral constraint on the total 
ionizing photon production of the Universe, this measure- 
ment is well suited to validate our models. Related to this, 
we briefly consider the contribution of metal-free stars to the 
cosmic near infrared background (CNIRB), as an excess in 
the range 1-2 /im has been discovered after subtraction of 
all known foreground sources (Cambresy et al. 2001; Santos, 
Bromm & Kamionkowski 2002). 

Another observational probe with great promise con- 
cerns the distribution of metals in the Universe, as con- 
siderable improvements have been made in measuring the 
metallicity of the IGM by quasar absorption lines (e.g. Pet- 
tini 2001). By computing the mass fraction of metals ejected 
into the IGM, we compare our results to metal abundances 
of the LyQ forest (e.g. Songaila 2001), thereby elucidating 
the enrichment properties of Pop III and Pop 11.5. Finally, 
we consider how metal-poor Galactic halo stars can be ac- 
comodated within our theoretical framework (for a recent 
review see Beers & Christlieb 2005). 



3.1 Ionizing Photon Production 

The most significant change related to the release of the 
WMAP 3-yr results concerns the optical depth to Thomson 
scattering, as the refined measurement revealed a reduction 
by nearly 50 per cent. At first sight, this might imply that 
early star formation was greatly suppressed, yet due to revi- 
sions in other cosmological parameters as well (i.e. the fluc- 
tuation power at small scales), high-z structure formation is 
rather shifted to lower redshifts by Az ~ 5, implying that 
primordial star formation is still relevant (e.g. Alvarez et 
al. 2006b). To investigate this connection more thoroughly, 
we must determine reionization histories and optical depths 
produced by our SFRs. 

A straightforward manner of quantifying the evolution 
of H II regions is to assume that each photon escaping a host 
galaxy ionizes one intergalactic hydrogen or helium atom. 
Supplementing this by a term encompassing recombination 
of electrons results in a simple formula for the evolution of 
the volume filling fraction of ionized regions (e.g. Wyithe & 
Loeb 2003, and references therein): 



dQion(z) 

dz 



1 dnion(z) ^, , 

= ; aBTlbC (z) 

nt dz 



xQl^{z){l + zf 



dt 

dz 



(21) 



where nt is the comoving density in baryons, nion(z) the 
comoving density in ionizing photons, as the case B re- 
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combination coefficient and C{z) the clumping factor. Tfiis 
metfiod is well established in modelling the evolution of H ii 
regions, and we refer to Wyithe & Loeb (2003) for details. 
Throughout this paper, we choose ob ~ 2.6x 10~^^ cm^ s~^, 
assuming that ionized regions are photoheated to ~ 10'* K 
(e.g. Osterbrock & Ferland 1989). To concentrate on the 
key effects, we make the following simplifying assumptions: 
(i) We neglect ionizing photons emitted by quasars, (ii) We 
do not distinguish between helium and hydrogen, implying 
that helium can only be ionized once, (iii) We use a sim- 
ple analytic fit of the form C{z) — 1 + 9[{1 + z) /7]^ , with 
C{z) = const. = 10 for z < 6, to describe the evolution of 
the clumping factor. This is a compromise between contem- 
porary approaches, and we refer to Haiman & Bryan (2006) 
for further discussion. 

The crucial quantity in equation (21) is nion(z), as it 
links our SFRs to the ionizing photon production. We estab- 
lish this connection by assuming that the number of ionizing 
photons emitted per redshift interval is proportional to the 
SFR. Explicitly, the conversion from newly formed mass in 
stars to the comoving density of ionizing photons is given 
by: 



1 dnion(z) 



rib 



dz 



_i_an 

Pm fib 



/csc'7ion**(z) 



(22) 



where /esc is the escape fraction of ionizing photons from 
the host halo and ?7ion the number of ionizing photons emit- 
ted per stellar baryon. Evidently, /esc and r;ion have different 
values for each population, and much effort has been under- 
taken to constrain these parameters. Specifically, Alvarez et 
al. (2006a) have found the escape fraction of ionizing pho- 
tons in minihaloes to be between 0.5 and 1, leading us to 
adopt /esc ~ 0.7 for Pop III. The escape fraction for typical 
Pop I/II galaxies is a strong function of redshift (e.g. Inoue, 
Iwata & Deharveng 2006), yet for early galaxies a reason- 
able estimate is 0.1, whereas for Pop II. 5 we can only guess 
that the best-fitting value must lie between 0.7 and 0.1. We 
hereby use /esc — 0.3, but keep in mind that it is generally 
a free parameter. In deriving T^ion for Pop III and Pop II. 5, 
we proceed similarly to Section 2, and use the spectra for 
metal-free stars provided by Schaerer (2002). This yields 
r?ion ~ 9 X 10'' for Pop III, and r/ion ~ 3 x W for Pop II.5. 
A common estimate for the ionizing photon production of 
Pop I/II is r]ion ~ 4 X 10^. 

With these quantities at hand, we compute the ion- 
ization histories Qion(z) implied by our SFRs. Panel (a) in 
Fig. 5 shows the ionization histories for the three distinct 
SFRs for Pop III determined in Section 2, i.e. early (dotted 
line), intermediate (dashed line) and late (dot-dashed line) 
reionization, combined with intermediate chemical feeback 
for Pop II. 5 (although variations in K„ hardly affect our 
results). By comparing Fig. 5 with Fig. 2, we conclude that 
intermediate reionization is by far the most likely and self- 
consistent model, as especially late reionization permits a 
genuinely high SFR (see Fig. 2), which induces a large con- 
tribution to the ionizing photon production, and ultimately 
contradicts the proposed reionization history. This further 
confirms our assumption that Pop III must be terminated 
fairly rapidly. 

In a final step, we determime the optical depth to 
Thomson scattering implied by these ionization histories: 
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Figure 5. Ionization histories (a), and implied optical depths 
(b) for the cases of preliminary early (dotted line), intermediate 
(dashed line) and late (dot-dashed line) reionization, as used to 
calculate the SFRs for Pop III (also see Fig. 2). For comparison, 
Panel (b) includes the range of allowed values for the optical depth 
as measured by WMAP (long-dashed lines). 



r = caTfii, 



dz' Qioniz') (l + z')" 



dt 



dz' 



(23) 



where z is the redshift of emission and ctt the cross-section 
to Thomson scattering (e.g. Wyithe & Loeb 2003). Accord- 
ing to this equation, the optical depth of a given ionization 
history is a well determined quantity, and Panel (b) in Fig. 5 
shows the optical depths corresponding to the three reioniza- 
tion histories of Panel (a). Not surprisingly, late reionization 
(as in Fig. 2) produces a value near the upper limit, whereas 
for early and intermediate reionization the resulting optical 
depth is very close to r = 0.09. 

At this point, we have also tested the validity of our star 
formation efficiencies for Pop III and Pop II. 5. For values 
much higher than our fiducial choices (see Table 1), we find 
optical depths that are unrealistically large. On the other 
hand, much smaller values completely neglect primordial 
star formation for ionization purposes. Taken together, this 
leads us to maintain the values we assumed in Section 2. 



3.2 Near Infrared Background 

The high ionizing photon yields of typical Pop III and 
Pop II. 5 stars might imply a significant contribution to the 
CNIRB. This possibility was investigated in detail by Santos 
et al. (2002), as an excess of the order ~ 1 nW m~^ sr~^ in 
the wavelength regime 1-2 /xm was detected by the Diffuse 
Infrared Background Experiment (DIRBE) on the Cosmic 
Background Explorer {COBE; e.g. Cambresy et al. 2001). 
The possible contribution to this wavelength range is made 
up of Lya photons, redshifted by the amount Az ~ 7 on the 
lower end, and Az ~ 15 on the higher end, implying that 
this period of time was relevant for producing the observed 
excess. Since Pop III is already extinct at these redshifts 
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(we only consider intermediate reionization) , Pop II. 5 is the 
only remaining candidate. To determine its importance for 
the CNIRB, we link the ionizing photon production to the 
mass in stars, and for simplicity assume that all ionizing 
photons are reprocessed through recombinations into Lya 



Jnir 



he 



ATTiriH 



r)\on 



*.W 



Az 



(24) 



Here, we have once again assumed that z/Nin/AfNiR is of the 
order of unity. Fairly independent of K^, this results in an 
upper limit of Jnir ~ 10~^ nW m~^ sr~^, which is much too 
low to account for the observed discrepancy. However, this 
result does not exclude the possibility of tracing primordial 
stars with fluctuations of the cosmic infrared background, 
as suggested by Kashlinsky et al. (2004). 
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Figure 6. Chemical Enrichment of the IGM by Pop III for inter- 
mediate reionization (dashed line), compared to the Lya forest 
and the critical metallicity Z^vit (long-dashed Unes). 



3.3 Metal Production 

In recent years, detailed studies of quasar absorption lines 
have provided a means of probing metal abundances of Lya 
systems with column densities as low as lO^** cm~^, made 
possible by improvements in data analysis techniques (e.g. 
Pettini 2001, and references therein). One of the most im- 
portant developments concerns the discovery of a metallicity 
floor in the Lya forest of the order 10~^ Zq in the redshift 
range 2; ~ 2 - 6, which cannot be explained by galactic out- 
flows of normal Pop I/II galaxies (e.g. Songaila 2001; Aguirre 
et al. 2005). This implies that early star formation might 
have been responsible for establishing such a homogeneous 
level of enrichment. Following up on this conclusion, we dis- 
cuss enrichment models for both Pop HI and Pop II. 5, and 
elucidate their contributions to the observed metal abun- 
dances of the IGM. 



3.3.1 Intergalactic Enrichment by Pop III 

As Heger et al. (2003) have discussed. Pop HI stars in the 
range 140 - 260 M0 end their lives as PISNe, extremely vio- 
lent SNe explosions with energies E^i < 100 and metal yields 
y < 0.5. Due to the shallow potential wells of minihaloes, 
these massive explosions completely evacuate their gaseous 
content, and distribute metals into the IGM with high ef- 
ficiency and homogeneity (e.g. Bromm et al. 2003; Yoshida 
et al. 2004). On the other hand. Pop HI stars with initial 
masses in the range 100 - 140 and 260 - 500 Mq collapse di- 
rectly to black holes, making them negligible for enrichment 
purposes (Heger et al. 2003). To find the correctly weighted 
metal yield, we average y over the IMF, resulting in y ~ 0.2. 
This quantity is essential in modelling the enrichment pro- 
cess by Pop HI, as it connects the mass in stars to the mass 
in metals. Applying the argument of homogeneous enrich- 
ment, we estimate a lower limit of the IGM metallicity by 
determining the ratio of density in metals to the background 
density for a given SFR: 

At 



Zigm{z) — — 

Pm 



*.(/) 



Az' 



Az' 



(25) 



Fig. 6 shows the resulting metallicity of the IGM for in- 
termediate reionization (dashed line) , compared to the fidu- 
cial value Z ~ 10~^ Zq for measurements of the Lya for- 
est (e.g. Songaila 2001), and the critical metallicity Zcrit = 



10" ■ Zq (long-dashed lines). For the case of highly effi- 
cient mixing, chemical enrichment by Pop HI does not raise 
the globally averaged metallicity above Zcrit, implying that 
Pop II. 5 star formation remains possible even in regions af- 
fected by PISNe winds. As is evident from Fig. 6, we pre- 
dict a minimum bedrock metallicity in the IGM of the order 
ZiGM ~ 10~^ Zq, much lower than current measurements 
of the Lya forest imply. This is not necessarily a contradic- 
tion, as the lowest column density regions of the IGM are 
not yet accessible for probing. Detecting such a low bedrock 
metallicity would provide a telltale signature of Pop HI star 
formation in early minihaloes. 



3.3.2 Local Enrichment by Pop II. 5 

The assumption of homogeneous metal distribution through- 
out the IGM breaks down for Pop II. 5, as typical explosion 
energies of conventional core-collapse SNe are significantly 
lower than those of PISNe, and the Universe at large has 
grown much more inhomogeneous since the epoch of Pop HI. 
We therefore assume that metals are distributed only inside 
enriching superwinds, and use equation (15) to determine 
the average density of such an enriched region: 



P(^) 



AMMri(M, z)nsT{M, z) 



(26) 



Inserting p instead of pni into equation (25), and adopt- 
ing the relevant values of y, E51 and /, for Pop II. 5, we 
determine the average metallicity of regions affected by su- 
perwinds. For our fiducial parameter K^ = 2/3 this yields 



2 X 10 Zq, which is in very 



agreement with 



our calculation of Section 2, and lies close to measurements 
of low column density Lya clouds (e.g. Songaila 2001). The 
metallicity exceeds Zcrit by almost one order of magnitude, 
further confirming our initial assumption that superwinds 
quench metal-free star formation in affected regions. An in- 
tuitively unexpected result is that stronger chemical feed- 
back results in less efficient metal production, which can be 
explained by the fact that a decrease in radiative losses en- 
ables the enriching winds to pollute larger regions of the 
IGM, thereby terminating Pop II. 5 more efficiently. 
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3.3.3 Metal-Poor Halo Stars 

The discovery of metal-poor stars in the halo of the Milky 
Way with extremely low iron content, yet excess light ele- 
ment abundances, has triggered discussion on the properties 
of the first stars (e.g. Beers & Christlieb 2005) . Various au- 
thors have argued that a less top-heavy IMF is more suc- 
cessfuU in reproducing the observed abundance patterns, as 
faint SNe and HNe of intermediate mass stars produce metal 
yields consistent with those found in metal-poor halo stars, 
provided these are of second-generation and formed from 
the relics of a metal-free population (e.g. Heger & Woosley 
2002; Umeda & Nomoto 2002, 2003, 2005; Tumlinson et al. 
2004; Daigne et al. 2004; Iwamoto et al. 2005). Based on the 
detailed physics of gas cooling, we have found primordial 
star formation to be dominated by stars > 10 Mq , and can 
therefore naturally account for the existence of metal-poor 
halo stars. 
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Figure 7. Three epochs of cosmological star formation history: 
Pop III (dotted line), Pop II. 5 (dashed line) and Pop I/II (dot- 
dashed line). 



4 SUMMARY AND CONCLUSIONS 

We have constructed star formation histories for two distinct 
populations of metal-free stars, characterized by their re- 
spective IMFs and formation environments. Specifically, we 
have attributed Pop III to form in minihaloes with charac- 
teristic stellar masses > 100 Mq , whereas second-generation 
haloes of > 10* M© give rise to an intermediate Popula- 
tion II. 5, enabling the formation of > 10 Mq stars by ac- 
tivating HD as an important low-temperature coolant. Ac- 
knowledging that both populations are subject to radiative 
and chemical feedback, we have shown that Pop III is ter- 
minated in the course of reionization, whereas Pop II. 5 ends 
with the distribution of metals into the IGM by enrich- 
ing superwinds, raising the metallicity in affected regions 
above Zcth = 10~^'^ Zq. We have confirmed the validity of 
our models by determining reionization histories and optical 
depths, and have excluded a dominant primordial contribu- 
tion to the excess of the CNIRB background in the wave- 
length regime 1-2 /im. Finally, we have discussed chemical 
enrichment of the IGM by Pop III and Pop II. 5, finding that 
Pop III produces a minimum bedrock metallicity in the IGM 
of the order Zjgm ~ 10~^ Zq, whereas Pop II. 5 enriches the 
IGM to levels of the low column density Lya forest. Summa- 
rizing these results. Fig. 7 compares the best-fitting SFRs 
for all three stellar populations, whereas Table 1 shows their 
characteristic quantities. 

The introduction of an intermediate generation of 
metal-free stars by means of formation environment pro- 
vides a possible physical explanation to the empirically de- 
rived argument that early star formation might have been 
dominated by a less top-heavy IMF (e.g. Umeda & Nomoto 
2002, 2003; Tumlinson et al. 2004; Tumlinson 2006a,b). One 
hint comes from the realization that a single, top-heavy pri- 
mordial population might overproduce metals while under- 
producing ionizing photons (e.g. Ricotti & Ostriker 2004; 
Daigne et al. 2004, 2005; Matteucci & Calura 2005). Some 
authors have also pointed out that metal yields of low col- 
umn density Lya clouds at 2; ~ 2 - 6 show enrichment sig- 
natures by more intermediate mass stars (e.g. Daigne et al. 
2004, 2005; Matteucci & Calura 2005). Another argument 
focuses on the abundance patterns of metal-poor Galactic 
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Table 1. Characteristic quantities of Pop III, Pop II. 5 and 
Pop I/II, from top to bottom: Typical stellar mass, minimum 
halo virial temperature, star formation efficiency, ionizing photon 
production, escape fraction, IMF-averaged explosion energy and 
IMF-averaged metal yield. 



halo stars, as they indicate prior enrichment by conventional 
core-collapse SNe instead of PISNe (e.g. Heger & Woosley 
2002; Umeda & Nomoto 2002, 2003, 2005; Tumlinson et al. 
2004). By additionally considering the chemical evolution 
of the Galaxy, it has been found that a more moderate 
IMF dominated by > 10 Mq stars best fits the currently 
available observational constraints (e.g. Tumlinson 2006a,b). 
Most recently, the revision of the optical depth to Thomson 
scattering has further weakened the case for a purely top- 
heavy primordial population, even in the light of a signif- 
icantly reduced amplitude of the matter fluctuation power 
(e.g. Choudhury & Ferrara 2006; Haiman & Bryan 2006). 
Taken together, these arguments provide a strong case for 
a more moderate IMF, yet thus far no physical framework 
has been found to accomodate such lower mass stars. In this 
paper, we have shown that the physics of HD cooling lead to 
a natural transition in terms of typical stellar masses, and 
force a distinction between two populations of metal-free 
stars. 

What was the relative importance of Pop HI and 
Pop II. 5? From our models, we estimate that only ~ 
10 per cent of all mass involved in metal-free star forma- 
tion went into Pop HI, and that Pop II. 5 was the dominant, 
and typical, mode of primordial star formation. This might 
suggest that Pop III had a negligible impact on observable 
quantities, yet due to their extreme nature they might have 
nevertheless left distinct signatures. Especially their influ- 
ence on the IGM could have been notable, as we expect 
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that only PISNe are able to enrich the lowest column den- 
sity regions in a homogeneous manner. Another clue might 
be provided by a characteristic reionization signature, as 
Pop III stars have unique spectra (e.g. Tumlinson & Shull 
2000; Tumlinson, Giroux & ShuU 2001; Bromm et al. 2001b; 
Oh, Haiman & Rees 2001; Tumlinson, Shull & Venkatesan 
2003; Venkatesan, Tumlinson & Shull 2003). In light of these 
findings, it is essential to push observational frontiers to ever 
higher redshifts, yet also extend the reach of numerical sim- 
ulations to more complex cosmological environments. In fu- 
ture work, we plan to apply such simulations to address the 
issues of intergalactic metal enrichment and star formation 
in atomic cooling haloes. 
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